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1 Introduction 


Early work of Axford and Hines ( 196 I) and Dungey (1961) on the configure* 
tlon of the earth's distant geomagnetic field due to the Interaction of the 
solar wind and the/^arth's geomagnetic field led to the now generally accepted 
view of tho g,ross features of the magnetosphere. The effect of the interac- 
tion is to generate current systems ^hich result in the earth's geomagnetic 
field on the sunward side being compressed and on the anti-sunward side being 
drawn out into a long tail-like structure. The magnetic field lines which 
form the tail of the magnetosphere originate from the polar regions of the 
earth's surface. One result of the extended tail of the magnetosphere la that 
there will be a continual plasma pressure gradient away from the earth paral- 
lel to the magnetic field lines between the relatively dense ionospheric plas- 
ma at about 300 km altitude and the distant geomagnetic tall. 

Initial stvidies of the effects of plasma density gradients on the diffu- 
sion dominated^ regions of the ionosphere had concentrated on applying plasma 
transport equations to solve for the case of the diffusive equilibrium distri- 
bution of mvjlti-ion species (cf- Mange, I960). In this case, the pressure 

'! 

gradient exiS'tlhg between the F-reglon peak and higher altitudes results in 
4he ions distributing thomselves such that the polarization electric field set 
up by the pressure gradient is just balanced by the gravitational force on the 
ions. Under these circumstances the outward flux of ions is zero. When the 
concept of open or :extended field lines became accepted, it was realized that 
the possibility existed for ions to continuously flow out from the earth into 
distant regions of the magnetosphere linked by geomagnetic field lines to the 
polar ionosphere* This condition is one of the solutions of the ion density 
distribution under diffusion dominated conditions and is known as 'dynamic 
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equilibrium* as opposed to 'diffusive equilibrium* discussed earlier. Both 
types of equilibrium represent steady state solutions to the plasma transport 
equations resulting in the requirement that the outward flux is a constant, 
zero in the case of diffusive equilibrium and some finite value in the case of 
dynamic equilibrium (Banks and Kockarts, 1973; Bauer, 1973)* Boundary condi- 
tions on the rate of production of the outflowing ions and Coulomb collisions 
with stationary ions result in an upper limit to the outward flux* The most 
dramatic ionospheric effect of dynamic equilibrium is to produce a marked re- 
duction in ion density of the outflowing species compared with the diffusive 
equilibrium situation. For a single ion species in its parent neutral gas, 
with electron, ion and neutral temperatures equal, the scale height is decre- 
ased by a factor of two for the change from diffusive to dynamic equilibrium. 

For a mixture of light and heavy ions such as and 0^, detailed solu- 
tions of the plasma transport equations show that the gravitational attraction 
on the heavier 0^ ion is such that the polarization electric field is not suf- 
ficient to permit the dynamic equilibrium condition and the ions form a diffu- 
sive equilibrium altitude-density profile. However, is light enough that 
once the inhibiting effect of collisions is reduced as the 0* density 
falls off, the polarization electric field can drive the ions into a continu- 
ous outflow parallel to the geomagnetic field (Banks and Holzer, 1968). This 
outflow of light ions from the polar regions of the earth's ionosphere was 
termed the polar wind by Axford. Further studies (Banks and Holzer, I969a» 
1 969 b) have shown that He^ is also light enough to flow out along open geomag- 
netic field lines, although the characteristics of the outflow differ consid- 
erably from that of H'*’, 
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Aa the polar wind Iona inoreaae their outflow velooity when they eeoape 
the effeota of Coulomb colliaiona with 0*, the conventional plaama tranaport 
equationa eventually become invalid due to the tranaltion from oolliaion domi- 
nated motion to oolliaionleaa motion* To accommodate the oolliaionleaa re- 
gion, kinetio models of the polar wind have been developed by Deaaler and 
Cloutier (1969) and by bemaire and Scherer (1970, 1971, 1972). Compariaona 
between the hydrodynamic and kinetic modela of the oolliaionleaa polar wind 
have shown that a good match occurs between the two techniquea (Holaer et al, 
1971). 

The studies of both the diffusive and dynamic equilibrium distribution of 
ion density with altitude which have been published are of steady state condi- 
tions. At present, no detailed studies of the time dependent nature of the 
outflow have been made, although some speculations on the characteristics of a 
shock front propagating along a suddenly depleted flux tube have been made 
(Banks et al, 1971). 

It is apparent that this outflow of plasma will provide a source of mag- 
netoBpherio plasma in the distant magnetospheric tail with its origin in the 
plasma generated in the polar regions by ionization of the earth's neutral at- 
mosphere by solar ultraviolet light in the sunlit pole and by particle ioniza- 
tion over the whole of the polar cap but concentrated in the auroral regions. 
The outflowing plaama into the magnetospheric tail is generated as a relative- 
ly cool thermal plasma. However, recent studies have shown that quite drastic 
modifications occur to the distribution functions of both the ions and elec- 
trons are predicted to occur as the plasma flows awhy from its source region 
in the ionosphere to great distances along the geomagnetic field lines. 

• If, 
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As was mentioned above, the ultimate aourqe of the polar wind ia the 
polar ionospheric plasma which has a peak density at around ?00 knt altitude. 
The source of to form the main element of polar wind is the accidentally 
resonant charge exchange reaction 

0"^ + H 0 

Thus, the amount of available to flow away from the earth is directly con- 
trolled by the ambient density of 0* ions in the topside ionosphere. Recent 
studies of the distribution of ion densities over the polar cap as a function 
of position and time hove shown that the 0^ density is by no means uniform. 
Hence, different flux tubes will have different initial plasma pressures to 
drive the polar wind, resulting in varying outward fluxes. The large varia- 
bility of 0'*’ density is a result of mapping the convection electric field 
about the magnetic pole which rotates diurnally about the geographic pole in 
the geographic reference frame (Sojka et al, 1979, 1981a, 1981b). 

The source of He^ ions flowing away from the earth' s ionosphere is less 
variable than the source, since it is produced primarily by chemical bal- 
ance between the photoionization of He and charge exchange between He“^ and the 

molecules 0^ and N 2 * However, local heating of the neutral atmosphere by 
plasma drift or by particle precipitation will increase the loss processes for 
He'*' in localized regions. The sources and Coulomb collision processes in the 
topside ionosphere for the H’’’ and He* ions result in a flux limited outflow 
having values as large 3 x 10®cm'^sec”^ for H* ions and 2 x lo'^cm'^sec”^ 
for He* ions. ' 

The fate of the outflowing plasma when it reaches the magnetosphere has 
been studied theoretically by Schunk A Watkins (1981, 1982) who have shown 
that marked anisotropies are predicted to develop as the polar wind plasma 
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flows away from th« earth* ITheae anlsotroplea may be unatable and lead to 
wave-partlole InBtabilltiea either from the flowing pl/|iisflia Itself or from an 

interaction with the environment through which the plasma it^ , flowing. 

. ( 

At present measurements of the low energy ion component of the magnetos- 
pherlc plasma have shown that unusual pitch angle distributions are not uncom- 
mon (Horowite et al» 1981) and that a significant population of ions with en- 
ergies around 10 eV and higher exists (Chappell, 1980). Unfortunately instru- 
mental problems associated With vehicle charging Often preclude ion measure- 
ments below 10 eV, and apart from an experiment flown on the ESA GEOS space- 
craft, there have been no Instruments flown which concentrate solely on 
measuring the detailed distribution function cf the low energy electron popu- 
lation in the magnetospheric plasma. 

In this study we had planned to use low energy ion distribution measure- 
ments from the ISEE spacecraft to guide the theoretical interpretation of 
processes occurring as the polar wind flows into the magnetosphere. However, 
we only received an nndocumented data tape near the end of the nominal period 
of performance. Conseciuently we started developing a new theoretical approach 
to predict the properties of low energy plasma in the magnetosphere f These 
studies were directed to studying the effects of wave-particle interactions 
involving the concept of plasmons and thereby drawing on quantum mechanical 
formulation for the processes occurring and also to study bulk energization of 
the low energy plasma through the concept of the energy momentum tensor for 
the plasma and its electromagnetic environment. 
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2. Comparison of Macroscopic and Microscopic Approach 

The dynamical behavior of high energy plasma may be determinad using sin- 
gle particle equations. For low energy plasmai collisional effects are impor- 
tant and one must use a fluid picture. This leads to a macroscopic view con- 
sidering mass, momentum and energy flows. These flows are related to macros- 
copic gradients by diffusion coefficients which are determined by microscopic 
processes such as wave- particle scattering. In addition, the effects of the 
macroscopic flows are related to measurements through distribution functions 
which are microscopic quantities. We, therefore, took a two - pronged appro- 
ach, both macroscopic and microscopic, to study the interaction of flowing 
plasma populations. 

2.2 Macroscopic Modeling Technique 

The conventional fluid equations are differential equations which are 
difficult to use with measured data because of noise in the measurements which 
can lead to erroneous derivatives. We use, instead, the general principle 
that a differential equation can be formulated as an integral equation (ince, 
1926*, Chandrasekhar, 1961). Thus, instead of working from the force equation, 
that is, the differential form of the momentum balance equation, we considered 
the momentum flux directly, that is, the integral form of the momentum balance 
equation. Further details of the technique are contained in a submitted paper 
entitled ''Energy and Momentum Flow in Electromagnetic Fie7<,ds and Plasma". 

By solving the integral form of the momentum balance equation, we deter- 
mined the momentum flow with several advantages over solutions of the differ- 
ential equation* 
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( 1 ) Our approach ia conceptualised differently and more ayatematlcally 
than the differential approach. We look at the energy-momentun tensor, spec- 
ifically at the components of this tensor which are the momentum fluxes. This 
Is conceptually different hecause now we realise that fields as well as parti* 
clea may carry momentum. The approach Is .more systematic because we can easi- 
ly handle momentum flow of fields and particles as well' as momentum flow due 
to friction. That Is, It is easier to identify the mechanisms which carry 
momentum and understand the contribution from each. 

(2) Differential equations are most commonly written in a given coordi- 
nate system for instance spherical coordinates. Our integral equations are 
written in a general coordinate system whose geometry is determined by the me- 
tric of the coordinate space. This along with General Relativistic notation 
greatly simplifies the algebra used to derive results. When our algebra is 
finished, we can apply our results in any coordinate system by choosing the 
metric for rectangular, spherical, cylindrical or any other coordinate space. 

(5) Data can be used without generating large relative errors inherent in 
numerical differentation. 

(4) The method is general and basic and can be used on almost any data. 
With proper statistical handling it can be applied to turbulent data. 

These four advantages to balancing the momentum are also present in the 
integral form of energy balance. Balancing the energy ia more difficult than 
momentum balance because we must use the thermodynamic characterisltics of the 
plasma which are not entirely known. 



f 


t 


Pag© 9 


Limited ua© ha© been made of this method by Yoh et. al» (1981) who der- 
ived a theoretical expreasion for the energy coupling of the aolar wind and 
the magnetoaphere. They looked at the apacial tranaport in the magnetic field 
by uaing part of the energy momentum tenaort Cowley (1980) hao made a atudy 
of raagnetoapheric energy tranaport where he conaidered the, balance between 
Joule heating and the apacial tranaport of energy by the electromagnetic 
field* Spacial energy flow in the magnetic fielda can flow into Joule heat- 
ing, but it can alao go into internal energy of the plaama or be uaed to ac- 
celerate the plaama. A large scale energy source such as low order multipoles 
generated by currents in the earth may provide a ’production range' energy 
source for MHD turbulence. Cowley (i960) writes down (an incomplete set) the 
momentum balance eq^uations but does not make use of them. 

We have formalized this work, Introducing a notation which is brief and 
convenient and to more fully utilize the power of the macroscopic approach. 

Even though most of the basic ideas are contained in classical electro- 
magnetic, Stratton (1941) or J.D. Jackson (l962), the moot complete diacua- 
sion, however, occurs in General Relativity, Weinberg (1972). General Rela-s 
tivity operates in 4*‘space with a given metric for this space resulting in a 
notation which is brief and manageable. By choosing the metric to be rectan- 
gular, spherical or cylindrical we can readily operate in any of these coordi- 
nate systems. 


The energy momentum tensor for a perfect fluid in a magnetic field is 
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+ u^u®) 4- pU®U® 

+ (F“ P®5 - 1/4 g«PpvMp ) 

5 v|i 

p ia the proper proaaure; here ia the ith component of proper 
velocity end ia ^ " 0 - j o,0 » (0,1, 2, 3); i - (1,2,3) and veloci- 
tiea have been normalized to the apeed of light. 

This ia the integral form of the momentum and energy balance equation, 
integrating T“P over 4 - apace producea a conatant. The differential form of 
the MHD equationa are contained in T®P by taking the 4 - divergence and equat- 
ing it to zero, T®B, • 0. 

0 

In the paper mentioned above we have been careful to take limita in the 
correct order, which ia not alwaya done,* Mianer et. al. (1975)* We took the 
Y ^ 1 limit after taking the derivatives. In this way we derived from this 
equation the hydrodynamic equations of Landau and Lifahitz (1959), and the 
magnetohydrodynamic equations of Lundquist (l952). 

At the present stage of development, we have formulated the situation for 
an ideal fluid in an electromagnetic field (Parish & Raitt, Energy and momen- 
tum Flow in Electromagnetic Fields and Plasma, Planetary and Space Sciencea , 
submitted, 1981a). The next step is to include a viscous tera in the energy 
momentum tensor allowing us to consider the momentum and energy flow due to 
viscosity in a dissipative fluid in an electromagnetic field. 

Satellite data requires special care in uae, beyond proper data reduction 
to remove instrument errors, even ideal measurements from a satellite are one 
point in space and time (ISEE A 4 6 make two point measurementa but the argu- 
ments are the same). If the velocity field, the spacecraft is measuring, var- 
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l9B Blovly in time and in space, then the neasurenant made represents the 
average velocity at that point. If, however, the velocity field is strongly 
turbulent, the velocity field is a random variable and the single measurement 
made could differ significantly from the actual average. 

In addition to applying our method to quiet areas, we have set up to make 
a proper statistical analysis of data in turbulent areas. This involves find- 
ing the correlation length and determining the correlation between velocity 
components, magnetic field components and between velocity and magnetic field 
components. Because of the bilinear nature of the terms in the energy momen- 
tum tensor which may be evaluated from two point correlation functions, and 
because the tensor deals with momentum and energy flows, it is necessary to 
investigate the relation of our new procedure to older turbulent transport 
theory of Kraichnan (I976a,b), Kraichnan and Nagrajan (1967) and Fyfe et. al. 
{ 1977 ) • Strong turbulence is non-linear and the large scales of planetary en- 
vironments make satellite measurements of turbulence unique and valuabli^. 
This application holds the possibility of testing and refining the MHI) turbu- 
lence theories. 

2.3 Microscopic Modeling Technique 

Following the observations of outwardly accelerated ions possibly gener- 
ated hy the mechanism which Inwardly accelerates auroral electrons, we worked 
on a model which studies the interaction of oppositely streaming plasmas to 
see if the heating observed by various workers (cf. Chappell (i960), Chappell 
et al. (I 98 O)) referred to earlier is predicted. Our model predicts the dis- 
tribution of the resulting plasma for direct comparison with the distributions 
measured in situ by magnetosphere satellites such as the ISBE-1 low energy 
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plasma analyzer or the BSA 6B0S experiment which also has low energy plasma 
analysis included in the payload. 

There are two commonly used methods of handling a plasma numerically when 
it is considered as a collection of individual particles rather than a fluid. 
Blther one follows the exact orbits of the particles resulting in the tracking 
of 3H coordinate points for N particles or one uses a statistical or Monte 
Carlo method and group the particles into velocity bins by means of a velocity 
distribution function and then uses the Vlasov equation. The former approach 
has the disadvantage of needing very large computer resources and also gener> 
ating so many instabilities that they are hard to characterize. The latter 
approach uses less computer resources, but tends to be too restrictive in the 
solutions attainable. The approach we have adopted falls between the two com- 
monly used methods discussed above. This technique allows a variety of insta- 
bilities and non-linearities to be handled in a controlled manner and the ef- 
fective wave-particle iuueraotions occurring to be quantified. 

At present, we have limited ourselves to Langmuir wave-particle interac- 
tions to improve our understanding of the model. Our Initial experience has 
indicated that this modeling approach will require less computer resources 
than other methods and will also be readily extendable to include other forme 
of plasma wave instability than the Langmuir waves we are currently studying. 

Wave-particle interactions dominate binary collisions in the magnetos- 
phere and the approach we are developing is a new method of studying space 

u 

plasmas and involves the application of Quantum Statistical Mechanics to these 
plasma interactions. The general basis for this technique has been described 
by Harris (1969), Pines (1962), Pines and Schrieffer (1962), Abrikosov, Gor- 
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kov, DryaloBhinaki (I963)p F«tt«r and ¥alaaka (l97l) and jpaftranoaa thareln. 

The phenomena we are inveatigating are claBaical with no quantum effecta, 
and it ifl poBBihle to develop thia formulation along olaflaioal linea (See Mar- 
tin, Siggia and Roae (1973); Krommea (1979); Boae (1979); Miahima, Petroeky, 
SuBUki (1979))« Thia claaaical formulation, however, ia in oarly atagea of 
development while the quantum mechanical theory ia highly developed and ao we 
take advantage of thia hy uaing the quantum mechanical treatment, in the claa- 
sical limit h ^ 0* Thia theory ia very powerful, and ia well auited to handle 
the large number of degreea of freedom and the nonlinearitiea inherent in 
wave-particle interactiona. 

Non - linearity in a plaama may be put into two categoriea. The firat ia 
an effective increaae in the maaa of the particles in the plaama due to the 
electromc'gnetic interactions of the charged particle with other particles. 
Correctins for this effective mass is know as *mass renormalization* • The 
second category is the correction of the electromagnetic force between two 
charges due to other charges being present; this is known as 'vertex renormal- 
ization* . Mass and vertex renormalization are not independent and plasmas are 
usually treated uaing 'vertex renormalization*. This allows us to account, 
for example, for the effects of a large amplitude wave on the plaama throxigh 
which it propagates, which in turn affects the wave itself. 

We start with the Schrodinger equation. Since the Vlasov equation may be 
derived from this equation (see von Roos (i960), Saenz (1957), Takabayasi 
(1954), Ross and Kirkwood (1957), Dahler (1959), Klimontovich (195S)), all so- 
lutions to the Vlasov equation are contained in our approach. The Schrodinger 
equation is first order in time and along with initial conditions will give 
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the time development of a quantum mechanical eyatem* The aye tern la deacribed 
by a wave function which givea the atatea of the particlea and wavea which 
make up the ayaten. 

In formulating the plaama proceeaea we aecond quantize » that la, we con- 
alder the electrostatic or electromagnetic field aa being made up of parti- 
cles. In the case of Langmuir waves , these particles are called plasmons. 

Ve represent the state of thr electrons or plasmons as 
t ■ J N(r)4j, 

where is the state of one particle, characterized by its momentum and ener- 
gy and N(r) is the occupation number or number of particles in the state a , 
We use Fermi's Golden Rule to find the change in the occupation number of a 
state r. The occupation number is the number of electrons or the number of 

plasmons in a state p , 

r 

3 N(r)/gt - 2ir/h 

( 1 ) 

2 

pg is the density of states around state s. |Mj,gj is the transition proba- 
bility of a transition from a state s to a state r. 

N can represent the niamber of electrons or the number of plasmons in 
state r. For discussion, we let N(r) be the number of electrons in state r, 
equations similar to those are deriving hold for plasmons. We now expand the 
right hand side of (l) into two terms. The first term represents the transi- 
tions which put an electron into state r and a second term which represents an 
electron in state r making a transition but of state r into any other state. 
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aN<£) ■ (tranaltiona per unit tlae Into atate r) 

3t 

- (tranaltiona per unit tine out of atate r) 

The number of tranaltiona per unit time Into atate r from another atate 
8» la the probability that the electrona will make the tranaltlon from a to r 
tlmea the probability that there la an electron In atate a tlmea the probabil” 
Ity that atate r la empty. State r muat be emptv for electrona alnce they are 
Fermlona and each atate cannot hold more than one. In our caae, an electuon 
nakea a tranaltlon Into or out of a atate by emitting or absorbing a plaamon. 

We may repreaent a single tranaltlon with a Feynmann diagram. For dla- 
cusaion, we consider only an electron In atate a emitting a plaamon Into a 
state t causing the electron to go Into state r. During the transition momen- 
tum Is conserved, and we use the momentum to label the states; state a la la- 
beled by p, state t Is labeled by k and state r la labeled by p-k. Time in- 
creases upward In this diagram: 



The point at which the transition takes place la called a vertex. We see 

2 2 

that associated with each vertex la the transition probability jM | . iM } 
la the probability that an electron emits a plasmon. The transition probabil- 
ities given by Harris (1969) or may be derived along the lines 

given by Bardeen and Pines (1955). The transition probabilities given by 
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Harris are for a bare vertex* We renomallse the vertex, which allows us to 
account for close Coulomb collisions. The long-range interactions, outside 

it 

the Debye sphere, are represented by the plasmone. If k is negative, the pro- 
cess Just described is the absorption of a plasmon by an electron. We now see 
that by generalizing this procedure to include all possible transltlcns, we 
can calculate the change in electron and plasmon distributions. We start with 
an initial distribution of electrons and plasmons and lot the computer keep 
track of all transitions and changes in the distribution functions. 

Up to now we have considered two beams streaming toward each other. This 
may, for instance, simulate two beams coming from conjugate hemispheres into 
the equatorial plane Just outside the plasmasphere. We consider only the in- 
teractions of the streaming electrons and each beam i@ assumed to have a neu- 
tralizing ion background with which it has no interactions. The two beams are 
assumed to be symmetric; that is, they have the same dfmsity, average flow 
velocity and temperature. 

We can obtain both the real and imaginary frequency; that is, the energy 
and decay rate of plasmons from this analysis. However, initially, we calcu- 
lated only the imaginary part of the frequency. The real part of the frequen- 
cy we find from standard linear analysis, singling out Langumir waves. This 
allows us to confine our attention to only one set of waves and look at their 
production and damping without the added complication of many different modes. 

The electron distribution was quantized into momentum spectrum steps of 
10*^ cm"^ and used 600 levels in the simulation. The simulation was one dimen- 
sional with 300 positive momentum steps and 300 negative momentum steps. In 
each momentijm bin, the number of electrons is N^(p). 
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We aterted nith two MaxweXliana atreaming toward each other with velocity 
oorreaponding to an energy of 10 eV, meaaured in units of m v^/2 relative to 

an inertial frame, ia the maaa of an electron. The temperature of each 
beam was 4 eV, and its density was 100 eleotrona/oo. The plaama frequency waa 

Uu ■ .1 MHz, and the time step waa u„/lO or 10”^ second. The initial diatri- 
bution ia shown in Figure 1. We ran our aimulation for 19 time atepa, and the 
final electron diatribution is shown in Figure 2. These figure show the two 
Maxwelliana overlapping and the quasitization of putting N electrons into a 
bin. 

As a check on the validity of our simulation, we continually calculated 
the total electron density which, electrons being Fermions, should remain con- 
stant. We see from Figure 3 that the total number of eloctrona waa conserved 
by the simulation at each time step. 

Since we treat the plasma waves as particles or plasmons, we have quan- 
tized their electric field. The plaemon distribution was quantized into mom- 
entum bins of size 10”^ cn“^ . The number of plasmons in bin p is Mx(p)» The 
simulation was started without plasmons and after 19 time steps, the plaaroon 
distribution developed to that shown in Figure 4* We see that eight bins, 
four with positive momenta and four negative bins, a reflection of the posi- 
tive ones through zero, contain plasmons. This shows the quantization of the 
waves as a certain number of plasmons contained in a bin. : 

The plasmons are bosons and their number is not conserved in the calcula- 
tion. Indeed, We have started with a plasma with no waves in it, and the sim- 
ulation has created Waves or plasmons. The time behavior of the total plasmon 
density is shown in Figure 5. ¥e see that the time behavior is very erratic 
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In the first few steps of generating the plasmons. 

These plasmons were generated hy spontaneous emission. As the number of 
plasmons increases with time, we expect stimulated emissions to become impor- 
tant. The simulation allows the plasmons to alter the distribution of the 
electrons. The altered distribution is used in the next timestep to generate 
plasmons. In this way the plasmons present in the plasma effect their own 
distribution. We never reached sufficient plasraon density for stimulated em- 
ission, in this case of only 19 timesteps, the plasmon levels were not large 
enough to alter the electron distribution (see Figures 1 and 2). 

In order to study the lifetime of the Langmuir waves corresponding to the 
decay of plasmons used in the model, we calculated the damping decrement given 
by 

y - 1/2N (p) gN (p)/gt 

XX 

at each timestep. The initial damping decrement is shown in Figure 6. There 
are no plasmons at the first timestep and so y was calculated from Harris 
(4.17) and corresponds to the Landau damping decrement. 

After 19 timestep, the damping decrement is shown in Figure We see 
that the plasmons generated are heavily damped aa shown by the 'spikes' great- 
ly exceeding the initial landau damping. 

We have shown that by using computer techniques on an old formalism that 
we can obtain reasonable results. We have produced waves from two counter 
streaming electron beams. We have followed the time development of these 
waves, and we have calculatedi the damping decrement, have prepared this new 
method to rederive the classical results of calculations and niimerical Simula- 
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ticn for the liraita on atability, wave apectrun and damping. 

We aet up the aoftware to continue the development of thia model and use 
a combination of the particle diatribution functiona and the wave fields meas- 
ured directly by magnetospheric satellites to test the predictions of the 
model and, thereby, establish the, validity of this new theoretical approach to 
modeling the dynamics of space plasmas. 

Thia method has been enormously successful in Quantum Electrodynamics in 
describing energy shifts (Lamb shift) and in Quantum Statistical Mechanics ap- 
plied to metals and super - conductivity. This is a very powerful method to 
handle the non - linearities which ariiie in particle - wave and wave - wave 
interactions. However, use of this technique in space plsama physics is in 
its infancy. The theory is readily applicable to a tenqus plasma, to a plasma 
containing a magnetic field, to a plasma containing large amplitude waves, to 
name just a few areas of research which will benefit from the research. 



5* Suggeationo for Future Vork 
5»1 Foundation for Future Work 
3 . 1.1 Maoroaoopic 

The work done on both the niiorosoopic and roacroacoplc theory have layed 
the foundation for future atudiea. 

The tnacroBoopic study is enormously powerful and particularly well-suited 
to study the large aoale measurements of Space Physics. This approach has 
been touched on but never used to its full extent. There seems to bo two be- 
liefs in Space Physics; firstly, that the macroscopic approach is unconven- 
tional and therefore inappropriate; and secondly, that all that is Involved is 
to use the energy momentum tensor. This last belief shows a lack of under- 
standing that there are ^several other aspects as discussed earlier. 

The macroscopic approach also leads us naturally into consideration of 
MHD turbulence. There is an enormous turbulence literature in the Journal of 
Fluid Mechanics and Physics of Fluids. These studies consider both hydrody- 
namical and raagnetqhydrodynamical turbulence and consist mostly of computer 
simulations. Turbulence depends on a parameter, the Reynolds number, which 
depends linearly on the size of the system. Turbulence is non-linear and 
Space Physics measurements are the only ones which have been made on a large 
enough system that internal turbulence measurements are available to test the 
computer simulations. 

The high time resolution magnetic field and particle distribution meas- 
urements now being made on a number of geophysical spacecraft are particularly 
well suited to a macroscopic study. With these we should be able to look at 
the turbulence in the magnetic and velocity fields and compare results to the 
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computer siraulatione. 

3.1.2 Microacppic 

Work on the macroscopic Is now suspended because of lack of funding; on 
the other hand, the microscopic theory which was designed to describe the 
polar wind is now being used to analyse the beam plasma experiments performed 
both in vacuum chambers and aboard the Space Shuttle. The microscopic theory 
applied to the magnetosphere requires wave data and particle data from current 
magnetospheric satellites, at present this work is also suspended due to lack 
of funding. However, we feel that the method holds promise for being able to 
model some of the wave particle interactions in a simpler manner and over a 
wider range of wave modes than is possible with quasi-linear or non-linear 
plasma theory* 
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6. PigurQ Captions 

Pigure t. Elactron distribution of oounter-atraaming eieotrona with 
at reaming 

velocity of 10 eV and temperature of 4 eV veraua momentum »t first 
timestep. 

Figure Electron distribution of oounter-oVrenming eleotrono with 
streaming velocity of 10 eV and temperature of 4 eV versus 
momentum at timestep 19» 

Figure 3. Total electron density versus time. 

Figure 4. Plasmon distribution versus momentum at timestep 19* 

Figure 5. Plasmon density versus time. 

Figure 6. Plasmon damping decrement versus momentum at first timestep. 
Figure ?• Plasmon damping decrement versus momentum at timestep 19* 
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